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Abstract 

We apply the Anharmonic Oscillator Symmetry Model to the description of vibrational excitations 
in and Td molecules. A systematic procedure can be used to establish the relation between the 
algebraic and configuration space formulations, by means of which new interactions are found in the 
algebraic model, leading to reliable spectroscopic predictions. We illustrate the method for the case 
of I?3 h -triatomic molecules and the T d Be-cluster. 
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The study of molecular vibrational spectra Q| requires theoretical models in order to analyze and 
interpret the measurements 0. These models range from simple parametrizations of the energy levels, 
such as the Dunham expansion to ab initio calculations, where solutions of the Schrodinger equation in 
different approximations are sought || |], || . In general, the latter involve the use of internal coordinates 
and the evaluation of force field constants associated to derivatives at the potential minima. While this 
method can be reliably applied to small molecules 0, it quickly becomes a formidable problem in the 
case of larger molecules, due to the size of their configuration spaces. New calculational tools to describe 
complex molecules are thus needed. 

In 1981 an algebraic approach was proposed to describe the roto-vibrational structure of diatomic 
molecules Q , subsequently extended to linear tri- and four- atomic molecules || and certain non- linear 
triatomic molecules Although these were encouraging results, the model could not be extended 
to polyatomic molecules, due to the impossibility of incorporating the underlying discrete symmetries. 
This difficulty could be surmounted by treating the vibrational degrees of freedom separately from the 
rotations. In 1984 Van Roosmalen et al proposed a U(2)-based model to describe the stretching vibrational 
modes in ABA molecules pi later extended to describe the stretching vibrations of polyatomic molecules 



such as octahedral and benzene like molecules 11 . Recently the bending modes have also been included 
in the framework, which was subsequently applied to describe C2u-triatomic molecules |l2| and the lower 
excitations of tetrahedral molecules (lj|, using a scheme which combines Lie-algebraic and point group 
methods. In a different approach, it has also been suggested to use a U(k + 1) model for the k — 3n — 3 
rotational and vibrational degrees of freedom of a n-atomic molecule. This model has the advantage that 
it incorporates all rotations and vibrations and takes into account the relevant point group symmetry 
10 , but for larger molecules the number of possible interactions and the size of the Hamiltonian matrices 
increase very rapidly, making it impractical to apply. 

Although the algebraic formulations have proved useful, several problems remained, most important 
of which is the absence of a clear connection to traditional methods. On the other hand, a related problem 
is the lack of a systematic procedure to construct all physically meaningful interactions in the algebraic 
space. In this paper we show that both these issues can be resolved by means of a general model for the 
analysis of molecular vibrational spectra, the Anharmonic Oscillator Symmetry Model (AOSM). In this 
approach it is possible to construct algebraic operators with well defined physical meaning, in particular 
interactions fundamental for the description of the degenerate modes present in systems exhibiting high 
degree of symmetry. The procedure to construct them takes full advantage of the discrete symmetry of 
the molecule and gives rise to all possible terms in a systematic fashion, providing a clear-cut connection 
between the algebraic scheme and the traditional analyses based on internal coordinates, which correspond 
to the harmonic limit of the model. 

As a test for this approach we apply the AOSM to the Be4 cluster and to three Pah-triatomic molecular 
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systems, namely H^, Be3 and Na^". Since small molecules can in general be well described by means 
of ab initio calculations |l5|, [T(| , we emphasize the basic purpose of this work. We have established an 
exact correspondence between configuration space and algebraic interactions by studying the harmonic 
limit of the U(2) algebra. This general procedure not only allows to derive the interactions in the AOSM 
from interactions in configuration space, but can also be applied to cases for which no configuration space 
interactions are available. The X^-triatomic molecules constitute the simplest systems where degenerate 
modes appear and where the new interactions in the model become significant. In the case of Be4, a 
direct comparison with ab initio calculations will be presented. The application of these techniques to 
more complex systems, such as the methane molecules, is presently under investigation p7| . 

The model is based on the isomorphism of the U(2) Lie algebra and the one dimensional Morse 
oscillator 

H = --- %+ D { e-^-2e-*) , (1) 

whose eigenstates £ can be associated with U(2) D SO{2) states . In order to see how this isomorphism 
comes about, consider the radial equation 

l(--T r T + r 4 + r2 )^ r )-( N + 1 )^) > ( 2 ) 

2 \ r dr dr r z J 

which corresponds to a two-dimensional harmonic oscillator (in units where h = fi = e = 1) associated 
to a U (2) symmetry algebra []l9t By carrying out the transformation 

„2 



r z = (N + l)e _p 



equation (2) transforms into 



d 2 



(e~ 2 P - 2e-P) 



0(p) = -m 2 0(p) , (3) 



dp 2 ' V 2 

which can be identified with (1) after defining x — pd and multiplying by h 2 /2pd 2 , provided that 



h 2 



D ^ 8^ + 1 > ' (4) 

£ ^ ■ (5) 

In the framework of the U (2) algebra, the operator N corresponds to the total number of bosons and is 
fixed by the potential shape according to (4), while m, the eigenvalue of the 5*0(2) generator J z , takes 
the values m — ±N/2, ±(JV — 2)/2, . . .. The Morse spectrum is reproduced twice and consequently for 
these applications the m-values must be restricted to be positive. In terms of the U(2) algebra, it is clear 
from (3-5) that the Morse Hamiltonian has the algebraic realization 

" - ^ - Aj * ■ < 6 > 
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In addition, the U(2) algebra includes the raising and lowering operators J + , J_, which connect different 
energy states in (3), while the angular momentum operator is given by J 2 = jN(N + 2), as can be readily 
shown. 

The Morse Hamiltonian (6) can be rewritten in the more convenient form 

H = AH M = ^[(j + J_ + j_J + )-N] , (7) 

where we have used the relation J 2 = J 2 — \ {J+ J- + J-J+) and added the constant term ^-j— in order 
to place the ground state at zero energy. The parameters N and A appearing in (7) are related with the 
usual harmonic and anharmonic constants uj e and x e ui e used in spectroscopy M . To obtain this relation 
it is convenient to introduce the quantum number 

N 

v= — -m, (8) 

which corresponds to the number of quanta in the oscillator [jl9]. This is seen by substituting (8) into 
(7). In terms of v, the corresponding energy expression takes the form 

E M =~(N+l/2) + A(N + 1)(v+ 1/2)- A(v + l/2) 2 , (9) 

from which wc immediately obtain 

co e = A(N + 1) , 
x e u, e = A . (10) 

Thus, in a diatomic molecule the parameters A and N can be determined by the spectroscopic constants 
w e and x e w e . 

We now consider the Ui(2) D SUi(2) D SOi(2) algebra, generated by the set {G;} = 
{Ni, J+.i, J-.ii Jo,j}j satisfying the commutation relations 

[Ja,i,J±,i] — ± J±,i j [J+,i,J-,i] = 2J ,i , [Ni,J^,i\ = 0, (11) 

with fx = ±,0. As mentioned before, for the symmetric irreducible representation [TV^, 0] of Ui(2) one can 
show that the Casimir operator is given by |l9| J? = Ni{Ni + 2)/4, from which follows the identification 
ji = Ni/2. The SOi(2) label is denoted by rrij. 

In the algebraic approach each relevant interatomic interaction is associated with a J7j(2) algebra 
As a specific example, we consider the Be4 cluster, which has a tetrahedral shape. V^h molecules can 
be similarly treated. In the Be4 case there are six Ui(2) algebras involved (i — 1, . . . , 6). The operators 
in the model are expressed in terms of the generators of these algebras, and the symmetry requirements 
of the tetrahedral group Td can be readily imposed 13 E(J. The local operators {G{\ acting on bond i 
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can be projected to any of the fundamental irreps T — Ai, E and F%. Using the J^.i generators ( |11| ) we 
obtain the Td tensors 

6 
i=l 

where fi = ±, and 7 denotes the component of T. The explicit expressions are given by 

6 

a/6; 



f A i - — V J ■ 



i=l 



1 

3 



F 2 



T% = ^(4,3-4,5) • (13) 

The Hamiltonian operator can be constructed by repeated couplings of these tensors to a total symmetry 
Ai, since it must commute with all operations in Td- This is accomplished by means of the 7rf-Clebsch- 
Gordan coefficients 13, ^0, 21 1. 

All calculations can be carried out in a symmetry-adapted basis, which is projected from the local 
basis 

t/i(2) ® ••• ® U 6 (2) D SOi{2) ® ••• ® S0 6 {2) D 50(2) 

I II II (14) 

I [Ni] , ... , [N 6 ] ; Vl , ... , v 6 ; V ) 

in which each anharmonic oscillator is well defined. By symmetry considerations, Ni — N for the six 
oscillators, i>i = Ni/2 — m,i denotes the phonon number in bond i and V — Vi is the total number of 



phonons |13|, 19 . The one-phonon states V — 1 are denoted by \i) with Vi = 1 and Vj^n = 0. Using the 
same projection technique as for the generators (|l3|), we find the six fundamental modes 



i=l 



£<*io- (15) 



The expansion coefficients are the same as in (p"3|). The higher phonon states v ' can also be constructed 
using the Clebsch-Gordan coefficients of Td |l^, [2(J . Since all operators are expressed in terms of powers 
of the Ui{2) generators, their matrix elements can be easily evaluated in closed form. The symmetry- 
adapted operators ([l3]) and states (|l5|) are the building blocks of the model. Note that for more complex 
molecules, the method allows the exact elimination of spurious states |17||. 
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We now proceed to expicitly construct the Be4 Hamiltonian. For interactions that are at most 
quadratic in the generators the procedure yields 

H = a\ Ha x + a,2 He + a 3 Hf 2 + b 2 Ve + b 3 Vf 2 , (16) 

with 



1 

2N 



= ( ^£ ' Tl + f 'I ■ f £ 



Vr = ^T r • f r , (17) 



Note that we have not included Va x in (|lq), since the combination 

£(Wr + Vr) = ~(tf + 2), (18) 
r 

is a constant. The five interaction terms in Eq. ( |l6| ) correspond to linear combinations of the ones 
obtained in lowest order in ]ll| , [l3] |. However, it is necessary to include interactions which are related to 
the vibrational angular momenta associated with the degenerate modes E and F% . These kind of terms is 
absent in the former versions of the model Wc 

now proceed to show how they can be obtained in 
the AOSM. In configuration space the vibrational angular momentum operator for the E mode is given 
by H 



^ 2 = -* Uf ^f-<12^e) , (19) 



L g B g 

where gf" and gjf are the normal coordinates associated to the E mode. This relation can be transformed 
to the algebraic space by means of the harmonic oscillator operators 



to obtain 

= -i(b B x h$-b^bf) . (21) 

Here b^ = J2i bi, with a similar form for w hile the can be read from (13). In order to find 
the algebraic expression for l A2 we first introduce a scale transformation in (11) 

fit = J^i/y/Ni , h = J+j/y/Wi . (22) 

The relevant commutator can be expressed as 

Ml - ^[W-.<] = ^2J , = 1-f , (23) 



where 



(24) 
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The other two commutators in (11) are not modified by (22). In the harmonic limit, which is defined 
by Ni — > oo , Eq. (23) reduces to the standard boson commutator [5j,5|] = 1. This limit corresponds 
to a contraction of SU(2) to the Weyl algebra and can be used to obtain a geometric interpretation of 
AOSM operators in terms of those in configuration space. In the opposite sense, Eq. ( p2] ) provides a 
procedure to construct the anharmonic representation of harmonic operators through the correspondence 
b\ — > b\ = J-.i/VN~i and b% — > b~i = J+^/yOVf. Applying this method to the vibrational angular 
momentum (|2l]) we find 

l M = -^(jV*2-^V£i) • ( 2 5) 

For the vibrational angular momentum l^ 1 associated with the F% mode we find a similar expression. 
The AOSM form of the corresponding Hamiltonian interactions is 

Hx = c 2 P 2 P 2 + c 3 *T l i" > (26) 

7 

With this method we can obtain an algebraic realization of arbitrary configuration space interactions. As 
a simple example, a one-dimensional harmonic oscillator Hamiltonian Hi = l/2(b\bi + bib\), transforms 
into 

i( J_,i J +ii + J+^J-.i) = i(J? - J 2 , ( ) = «, + 1/2 - ^ , (27) 

where in the last step we used relation (24). The spectrum of (27) has an anharmonic correction, analogous 
to the quadratic term in the Morse potential spectrum. We are thus substituting harmonic oscillators by 
Morse oscillators in the AOSM. 

A more interesting application is to use our model to fit the spectroscopic data of several polyatomic 
molecules. In the case of Be4 the energy spectrum was analyzed by ab initio methods in Jig ], where 
force-field constants corresponding to an expansion of the potential up to fourth order in the normal 
coordinates and momenta were evaluated. We have generated the ab initio spectrum up to three phonons 
using the analysis in ]2^] . For the algebraic Hamiltonian we take 

H = H + c 3 J2 tf 1 Z T + d 3 {Uf 2 + Vf 2 ) 2 + e 3 6 33 , (28) 

7 

where Hq was defined in (16). The term 33 is the algebraic form of the corresponding i^-mode in- 



teraction in 1 22 which is responsible for the splitting of the vibrational levels with the same I in the 
i*2 overtones. For T> 3 h molecules we can follow an analogous procedure, namely, we can construct the 
T> 3 h symmetry-adapted operators and states corresponding to (13) and (15) and carry out the building 
up procedure to construct the Hamiltonian and higher phonon states, using in this case the appropriate 
projection operators and Clebsch-Gordan coefficients |l^, |2Cj. Here we omit the details for lack of space 
and only present the fit to the energy spectrum (|^, [24j . 
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Note that the Be4 Hamiltonian (28) preserves the total phonon-number V. This is a good approxi- 



mation for this case according to the analysis of [122, 25 , but it is known that Fermi resonances can occur 
for certain molecules when the fundamental mode frequencies are such that (V, V) states with V ^ V 
are close in energy. These interactions can be introduced in the Hamiltonian but the size of the energy 
matrices grows very rapidly, so the best way to deal with this problem is through perturbation theory. 
We now present the results of our least-square fits to the energy spectra of Be4, Be3, Najj" and H^. 

In Table I we show the fit to Be4 using the Hamiltonian (28). The r.m.s. deviation obtained is 2.3 
cnr 1 , which can be considered of spectroscopic quality. We point out that in |^2|, several higher 
order interactions are present which we have neglected. Since our model can be put into a one to one 
correspondence with the configuration space calculations, it is in fact possible to improve the accuracy 
of the fit considerably, but we have lused a simpler Hamiltonian than the one of |22| ^5). When no 
ab initio calculations are available (or feasible) the AOSM approach can be used empirically, achieving 
increasingly good fits by the inclusion of higher order interactions [[l7]. In Table II we presenlt AOSM 
fits to the spectra of Be3, Na^ and H^" up to three phonons. While remarkably accurate descriptions of 
the first two molecules can be achieved using a four-parameter Hamiltonian, we were forced to include 
four additional higher order terms in the Hjj" Hamiltonian in order to properly describe this molecule. 
This is in accordance with the work of Carter and Meyer Jl6| , who were forced to include twice as many 
terms in the potential energy surface lfor H^" than for the Najj" molecule. The Hjj~ ion is a very "soft" 
molecule which, due to the light mass of its atomic constituents carries out large amplitude oscillations 
from its equilibrium positions |l6| ]. 

A finer test for the model is to use the wave functions to evaluate infrared and Raman transitions. The 
algebraic realization of the transition operators can be obtained from their expression in configuration 
space using the large N connection, or purely algebraically by their tensorial properties under the relevant 
point group f2~ij| . We remark that the model can also be extended to include the rotational degrees of 
freedom, by coupling the vibrational wave functions to rotational states properly symmetrized to carry 



the point group representations 22 



The AOSM is based on symmetry methods which systematically incorporate group theoretical tech- 
niques, providing a clear methodological procedure that can be applied to more complex molecules. We 
define symmetry adapted operators that have a specific action over the function space. This is a general 
procedure which gives rise to a clear physical interpretation of the interactions and has the additional 
advantage of considerably improving the convergence during the least square energy fits. Based on the 
harmonic limit of the SU(2) algebra we have found a systematic approach to derive an algebraic realiza- 
tion of interactions given in configuration space. The model surmounts one of the main objections raised 
against the use of algebraic models, where it was not possible to obtain a direct correspondence with 
the configuration-space approaches. For the general case when there is no information about the form 



of these interactions in configuration space, we have devised an algebraic procedure to derive them using 
their tensorial structure under the point group. The combination of the different methodologies leads to 
the AOSM, which can be applied in the same fashion to more complex molecules. 

We believe that the AOSM represents a systematic, simple but accurate alternative to the traditional 
methods, particularly for polyatomic molecules, where the integro-differential approaches are too complex 
to be applied or require very large numerical calculations. Since the model provides manageable wave 
functions, it is possible to evaluate the matrix elements of arbitrary physical operators, which have a 
simple representation in the algebraic space. The study of electromagnetic intensities, as well as the 
application of the model to more complex molecules is currently under investigation []l7], Q . 

This work was supported in part by the European Community under contract nr. CIl*-CT94-0072, 
DGAPA-UNAM under project IN105194, CONACyT-Mexico under project 400340-5-3401E and Spanish 
DGCYT under project PB92-0663. 
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Table I: Fit to ab initio calculations for Be4. The values of the parameters are N — 40, a\ — 646.95, 
a 2 = 463.51, a 3 = 675.54, b 2 = 129.46, 63 = 317.01, c 3 = -14.285, d 3 = 10.688, and e 3 = -1.948. The parameters 
and energies are given in cm -1 . (1/1, v™, v\) denotes the number of phonons in the Ax, E and F2, respectively, 
and m and I the value of the vibrational angular momenta pM. 
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Table II: Least square energy fit for H+, Be 3 and NaJ using the Hamiltonian (7.1). We show the energy differences 
AE = Eth — E exp . The values of the energies E exp are given in Table V. 
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